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This  paper  deals  with  finding  the  optimum  tilt  angle  of  solar  panels  for  solar  energy  applications.  The 
optimization  of  tilt  angles  was  performed  using  solar  radiation  data  measured  for  eight  big  provinces  in 
Turkey.  The  optimum  angle  for  tilted  surfaces  varying  from  0  to  90  in  steps  of  1°  was  calculated  by 
searching  for  the  values  of  which  the  daily  total  solar  radiation  was  at  a  maximum  for  a  specific  period. 
It  was  found  that  the  optimum  tilt  angle  changed  between  0  and  65  throughout  the  year  in  Turkey.  It 
was  seen  that  the  optimum  tilt  angle  reached  a  minimum  of  0  in  June  and  July  and,  the  monthly 
average  daily  total  radiation  at  this  angle  was  generally  at  a  maximum.  In  addition,  the  optimum  tilt 
angle  increased  during  the  winter  months  and  reached  a  maximum  in  December  in  all  provinces. 
Likewise,  general  correlations  were  developed  to  estimate  the  optimum  tilt  angle  of  solar  collectors 
used  in  Turkey  and  their  accuracies  were  compared  on  the  basis  of  statistical  error  tests  of  Mean  Bias 
Error  (MBE),  Root  Mean  Square  Error  (RMSE),  t-statistic  (t-stat)  and  correlation  coefficient  (r). 
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1.  Introduction 

Turkey  is  an  energy  importing  nation  with  more  than  half  of 
our  energy  requirements  met  by  imported  fuels.  Air  pollution  is 
becoming  a  significant  environmental  concern  in  the  country. 
Achieving  solutions  to  the  environmental  problems  that  we  face 
today  requires  long-term  potential  actions  for  sustainable  devel¬ 
opment.  In  this  regard,  renewable  energy  resources  appear  to  be 
the  one  of  the  most  efficient  and  effective  solutions.  For  govern¬ 
ments  or  societies  to  attain  sustainable  development,  a  great  deal 
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of  effort  must  be  devoted  to  utilizing  sustainable  energy  resources 
in  terms  of  renewable  energies  [1,2].  The  studies  on  the  new  and 
renewable  energy  sources  have  gained  speed  and  are  encouraged 
due  to  the  fact  that  energy  resources  used  today  run  out  rapidly 
and  cause  environmental  pollution  [3]. 

Estimates  of  the  monthly  average  solar  radiation  incident  on 
the  surface  of  various  orientations  are  required  for  solar  energy 
design  procedures.  Monthly  averages  of  the  daily  solar  radiation 
incident  upon  a  horizontal  surface  are  available  in  the  literature 
for  many  locations  [4-9].  However,  radiation  data  on  tilted 
surfaces  is  not  generally  available.  Tiris  et  al.  [5]  calculated  the 
correlations  of  the  monthly  average  daily  global,  diffuse  and  beam 
radiations  with  hours  of  bright  sunshine  in  Gebze,  Turkey.  Bakirci 
[7]  developed  the  correlations  for  estimation  of  daily  global  solar 
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Nomenclature 

K, 

Clearness  index 

H 

Monthly  average  daily  global  radiation  on  horizontal 

Rb 

Ratio  of  average  beam  radiation  on  a  tilted  surface  to 
that  on  horizontal  surface 

surface  (MJm~2-day) 

s 

Tilt  of  surface  from  horizontal 

Ha 

Monthly  average  daily  diffuse  radiation  on  horizontal 

$opt 

Optimum  tilt  angle  (°) 

surface  (MJm~2-day) 

S 

Solar  declination  (°) 

H0 

Monthly  average  daily  extraterrestrial  radiation 

cp 

Latitude  of  the  place  (°) 

(MJm~2-day) 

P 

Ground  reflectance 

HT 

Monthly  average  daily  global  radiation  on  tilted  sur¬ 

CDs 

Hour  angle  (°) 

face  (MJm~2-day) 

CD's 

Sunset  hour  angle  (°) 

radiation  with  hours  of  bright  sunshine  in  Turkey.  Zuhairy  and 
Sayigh  [8]  carried  out  simulation  and  modeling  of  solar  radiation 
in  Saudi  Arabia.  Ulgen  and  Hepbasli  [9]  investigated  the  diffuse 
fraction  of  daily  and  monthly  global  radiation  for  Izmir,  Turkey. 

The  amount  of  energy  collected  by  a  solar  collector  and  thus 
converted  to  thermal  energy  is  highly  influenced  by  its  orientation 
and  tilt  angle.  Designers  of  solar  systems  that  use  flat-plate 
collectors  need  to  know  the  amount  of  solar  energy  collected  by  a 
tilted  surface.  For  such  a  system,  the  collector  slope  is  one  of  many 
design  considerations.  Calculations  of  the  solar  energy  absorbed  by  a 
tilted  surface  directly  facing  the  equator  with  an  azimuth  are 
available  in  the  literature.  However,  designers  often  need  to  know 
how  much  energy  can  be  absorbed  by  the  surfaces  with  other 
azimuths.  It  is  generally  known  that  in  the  northern  hemisphere,  the 
optimum  collector  orientation  is  south-facing,  and  the  optimum  tilt 
depends  upon  the  latitude  and  the  day  of  the  year  [10]. 

Data  of  the  monthly  global  and  diffuse  solar  radiation  obtained 
from  observations  over  many  years  are  required  for  acceptable 
estimation  of  the  optimal  tilt  angle  of  a  solar  collector.  However, 
the  monthly  diffuse  radiation  is  not  always  available  in  many 
locations  in  Turkey.  Therefore,  in  this  article  an  empirical  correla¬ 
tion  suggested  by  Page  [11]  was  used  for  the  calculation  of 
monthly  diffuse  radiation  in  the  provinces  considered  and  a 
simple  mathematical  procedure  was  also  applied  to  estimate 
the  optimal  tilt  angle  of  a  solar  collector. 

The  main  objective  of  this  study  is  to  determine  the  optimal 
tilt  angle  of  solar  panels  used  in  Turkey  for  maximizing  their 
energy  collection  and  to  obtain  general  correlations  for  estimating 
the  optimum  slopes  that  are  precise  and  convenient  for  use. 


2.  Background  and  previous  works 

An  important  parameter  in  the  optimum  utilization  of  solar  panels 
is  their  tilt  angle  with  the  horizontal.  This  changes  the  amount  of 
solar  radiation  that  reaches  the  solar  panels.  The  optimal  tilt  angle  of 
a  solar  collector  is  related  to  the  local  climatic  condition,  the 
geographic  latitude  and  the  period  of  its  use.  Therefore,  different 
locations  will  have  different  optimal  tilt  angles  for  solar  collectors. 

A  number  of  studies  on  the  optimal  tilt  angle  of  solar  panels 
around  the  world  have  been  presented  [12-36].  However,  studies 
carried  out  in  Turkey  are  few  in  number  and,  the  general  correla¬ 
tions  giving  optimum  slopes  for  Turkey  are  not  obtained.  Iqbal  [12] 
investigated  the  optimum  collector  slope  for  residential  heating  in 
adverse  climates.  Soulayman  [15]  carried  out  a  study  on  the 
optimum  tilt  of  solar  absorber  plates.  Hussein  et  al.  [21]  developed 
a  theoretical  analysis  of  the  instantaneous,  daily,  and  yearly 
enhancement  in  solar  energy  collection  of  a  tilted  flat-plate  solar 
collector  augmented  by  a  plane  reflector.  Bari  [22]  described  a 
method  to  determine  the  optimum  slope  angle  and  orientation  of 
solar  collectors  for  different  periods  of  operation  at  any  position  in 
the  Malaysian  territory.  The  method  utilizes  both  the  direct  and 
diffuse  components  of  solar  radiation.  Tang  and  Wu  [25]  carried  out 


a  reasonable  estimation  of  the  optimal  tilt  angle  of  a  fixed  collector 
for  maximizing  its  energy  collection  based  on  the  monthly  global 
and  diffuse  radiation  on  a  horizontal  surface.  Gunerhan  and  Hepbasli 
[29]  calculated  the  optimum  tilt  angles  by  searching  for  the  values 
for  which  the  total  radiation  on  the  collector  surface  is  at  a 
maximum  for  a  particular  day  or  a  specific  period. 

Mehleri  et  al.  [31]  carried  out  a  study  on  the  determination  of  the 
optimum  tilt  angle  and  orientation  for  solar  photovoltaic  arrays  in 
order  to  maximize  the  incident  of  solar  irradiance  exposed  on  the 
array,  for  a  specific  period  of  time.  Moghadam  et  al.  [32]  performed 
optimization  of  solar  flat  collector  inclination.  Monthly,  seasonal, 
semi-annual  and  annual  optimum  tilt  angles  were  determined.  Ghosh 
et  al.  [33]  determined  the  seasonal  optimum  tilt  angles,  solar 
radiations  on  variously  oriented,  single  and  double  axis  tracking 
surfaces  at  Dhaka.  Three  mathematical  models  for  the  point  source 
with  parameters  optimized  for  a  variety  of  climatic  conditions  were 
employed  to  determine  hourly  and  seasonal  optimum  tilt  angles. 
Maatallah  et  al.  [34]  presented  an  overview  on  research  works  on 
solar  radiation  basics  and  photovoltaic  generation.  The  effects  of 
azimuth  and  tilt  angles  on  the  output  power  of  a  photovoltaic  module 
were  investigated.  Kaldellis  and  Zafirakis  [35]  carded  out  an  experi¬ 
mental  study  in  the  area  of  Athens  in  order  to  evaluate  the 
performance  of  different  PV  panel  tilt  angles  during  the  summer 
period.  The  angle  of  15°  (  +  2.5°)  was  designated  as  optimum  for 
almost  the  entire  summer  period.  Benghanem  [36]  performed  a  study 
on  the  optimum  slope  and  orientation  of  a  surface  receiving  a 
maximum  solar  radiation.  The  annual  optimum  tilt  angle  was  found 
to  be  approximately  equal  to  the  latitude  of  the  location.  Sirald  and 
Pillay  [37]  proposed  a  simple  method  on  a  modified  sky  model  to 
calculate  the  optimum  angle  of  installation  for  urban  applications.  It 
was  expressed  that  the  results  demonstrated  the  dependency  of  the 
optimum  angle  of  installation  on  the  latitude,  weather  conditions  and 
surroundings. 

Lave  and  Kleissl  [38]  calculated  the  optimum  tilt  and  azimuth 
angles  of  solar  panels  for  a  grid  of  0.1°  by  0.1°  National  Solar  Radiation 
Database  cells  covering  the  continental  United  States.  The  yearly 
global  irradiation  incident  on  a  panel  at  this  optimum  orientation  was 
compared  to  the  solar  radiation  received  by  a  flat  horizontal  panel 
and  a  2-axis  tracking  panel.  Lubitz  [39]  investigated  the  effect  of 
manual  tilt  adjustments  on  incident  irradiance  on  fixed  and  tracking 
solar  panels.  The  optimum  tilt  angle  for  an  azimuth  tracking  panel 
was  found  to  be  on  average  19°  closer  to  the  vertical  than  the 
optimum  tilt  angle  for  a  fixed,  south-facing  panel  at  the  same  site. 
Chau  [40]  performed  a  study  on  optimum  tilt  angles  for  solar 
collectors  in  clear  sky  conditions.  The  optimum  tilt  angles  for 
collectors  at  various  latitudes  were  determined  for  all  12  months  of 
the  year.  It  was  expressed  that  these  results  were  also  applicable  to 
collectors  with  a  cylindrically  curved  transparent  cover.  Moon  et  al. 
[41]  determined  the  optimum  tilt  angles  for  solar  collectors  in  171 
locations  in  North  America.  The  equations  that  predicted  the  opti¬ 
mum  collector  tilt  angle  were  derived  by  using  regression  analysis 
based  on  the  results  obtained  for  the  171  locations  for  several 
different  periods  of  the  year. 
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3.  Estimation  of  solar  radiation  on  an  inclined  surface 


The  monthly  average  values  of  solar  radiation  incident  on 
surfaces  of  various  orientations  are  required  for  solar  energy 
applications.  The  monthly  averages  of  the  daily  solar  radiation 
incident  upon  a  horizontal  surface  are  available  for  many  loca¬ 
tions.  However,  radiation  data  on  tilted  surfaces  are  generally  not 
available.  A  simple  method  to  estimate  the  average  daily  radiation 
for  each  calendar  month  on  surfaces  facing  directly  towards  the 
equator  has  been  developed  by  Liu  and  Jordan  [42]. 

The  monthly  average  daily  radiation  on  a  horizontal  surface 
( H ),  for  each  calendar  month  can  be  expressed  by  defining  KT,  the 
fraction  of  the  mean  daily  extraterrestrial  radiation  (H0) 

Kt  =  H/H0  (1) 


H  -  24l 
Ido  — 

71 


1  +0.033cos 


360n 

/365 


^cos<pcos<5sin<Ms  +  simpsind) 


(2) 

where  lsc  is  the  solar  constant  ( =  1353  W/m2),  n  is  the  day  of  the  year 
given  for  each  month,  q>  is  the  latitude,  3  is  the  solar  declination  and 
ojs  is  the  hour  angle.  Values  of  3  and  cos  are  given  by  [43,44]. 

'360(n  +  284f 


3  =  23.45sin 


(3) 

(4) 


365 

cos  =  arcos  [ — tan(<5  )tan(  qo )] 

The  average  daily  radiation  on  a  tilted  surface  (HT)  can  be 
expressed  by 

Ht  =  RH  (5) 


where  R  is  defined  as  the  ratio  of  the  daily  average  radiation  on  a 
tilted  surface  to  that  on  a  horizontal  surface  for  each  month.  R  can 
be  estimated  by  considering  individually  the  beam,  diffuse  and 
reflected  components  of  the  radiation  incidence  on  the  tilted  sur¬ 
face.  Assuming  the  diffuse  and  reflected  radiation  to  be  isotropic,  Liu 
and  Jordan  have  proposed  that  R  can  be  expressed  as  [42] 


R  = 


/l+cos(s)\ 
V  2  H  ) 


+p 


l-cos(s)^ 

2  J 


(6) 


where  Rb  is  the  ratio  of  the  average  beam  radiation  on  the  tilted 
surface  to  that  on  a  horizontal  surface  for  each  month,  s  the  tilt  of  the 
surface  from  horizontal,  and  p  (  =  0.2)  the  ground  reflectance.  Hd  is 
the  monthly  average  daily  diffuse  radiation  and  is  given  by  Page  [11] 

Hd  =  H(1.00-1.13KT)  (7) 


It  is  noted  that  Rb  is  a  function  of  the  transmittance  of  the 
atmosphere  which  depends  upon  the  atmospheric  cloudiness,  water 
vapor  and  particulate  concentration.  However,  Liu  and  Jordan  suggest 
that  Rb  can  be  estimated  to  be  the  ratio  of  extraterrestrial  radiation  on 
the  tilted  surface  to  that  on  a  horizontal  surface  for  the  month  [45]. 
For  surfaces  facing  directly  towards  the  equator 

^  cos((p-s)cos(<5)sin(co')  +  a)'  (71  / 1 80)  sin(qt>  — s)sin(<5) 

cos(<p)cos(d)sin(cus)  +  ct)s  (7t/180)  sin((p)sin(  d) 

where  mj  is  the  sunset  hour  angle  for  the  tilted  surface  and  is  given 

W 

c o's  =  min{a>s, arcos  [-tan(<p— s)tan(5)] }  (9) 


4.  Results  and  discussion 

The  present  study  was  carried  out  for  the  eight  main  provinces  of 
Turkey.  The  above  formulae  were  used  to  calculate  the  monthly 
average  daily  total  radiation  on  a  south  facing  surface  as  the  tilt  angle 
was  changed  from  0°  to  90’  in  steps  of  1°.  The  optimum  angle  for 
tilted  surface  was  calculated  by  searching  for  the  values  of  which  the 


daily  total  solar  radiation  was  at  a  maximum  for  a  specific  period. 
Therefore,  the  optimum  tilt  angle  was  found  for  certain  provinces  of 
Turkey. 

Information  for  the  provinces  considered  in  this  study  is  given  in 
Table  1.  The  recommended  average  day  and  declination  for  each 
month  are  presented  in  Table  2  [46],  Table  3  shows  the  measured 
monthly  average  daily  global  solar  radiation  for  these  provinces  of 
Turkey.  It  also  shows  the  calculated  diffuse  and  extraterrestrial  solar 
radiation  values,  clearness  index,  optimum  tilt-angle  and  monthly 
average  daily  global  radiation  on  optimum  tilted  surface  in  eight 
provinces  of  Turkey.  It  was  found  that  the  optimum  tilt  angle  changed 
between  0°  (June  and  July)  and  60’  (December)  for  Adana,  0°  (June) 
and  60’  (December)  for  Ankara,  0°  (June  ad  July)  and  61°  (December) 
for  Diyarbakir,  0  (June)  and  61°  (December)  for  Izmir  and  Trabzon,  0° 
(June)  and  65°  (December)  for  Erzurum,  0C  Gone)  and  59’  (December) 
for  Istanbul  and  0°  Quire)  and  62’  (December)  for  Samsun  throughout 
a  year.  The  monthly  average  daily  solar  radiation  availabilities  of 
tilted  surfaces  for  the  provinces  of  Turkey  are  presented  in  Fig.  1. 
Fig.  2  shows  monthly  average  daily  solar  radiations  for  the  deter¬ 
mined  tilt  angles  when  the  solar  collector  is  tilted  at  the  optimum 
angle,  the  seasonal  average,  the  yearly  average,  the  latitude  and  the 
latitude  ±15.  General  correlations  developed  for  the  optimum  tilt 
angles  of  the  solar  collectors  used  in  Turkey  are  given  in  Fig.  3.  Also, 
the  results  of  this  study  are  compared  with  the  studies  carried  out  by 
Ulgen  [28]  and  Gunerhan  and  Hepbasli  [29]  in  the  city  of  Izmir  in 
Fig.  4.  As  can  be  seen  from  Fig.  4,  the  results  obtained  from  each  study 
are  in  general  agreement. 

The  values  of  the  monthly  average  daily  global  solar  radiation 
for  the  provinces  of  Turkey  were  taken  from  the  Turkish  State 
Meteorological  Service  for  the  period  between  1991  and  2005. 
The  optimum  tilt  angles  were  determined  based  on  the  experi¬ 
mental  data  of  the  monthly  average  global  solar  radiation. 

Additionally,  the  correlation  equations  of  the  optimum  tilt 
angle  (Sopt)  based  on  the  declination  factor  are  developed  for 
Turkey  in  general.  These  correlation  equations  are  given  below 

Sopt  =  32.619-1.3629(5)  (10) 


Table  1 

Information  for  the  provinces  considered  in  the  study. 


Location 

Longitude  (£) 

Latitude  (N) 

Elevation  (m) 

Adana 

35.18 

36.59 

20 

Ankara 

32.53 

39.57 

894 

Diyarbakir 

40.12 

37.55 

660 

Erzurum 

41.16 

39.55 

1869 

istanbul 

29.05 

40.58 

39 

Izmir 

27.10 

38.24 

25 

Samsun 

36.20 

41.17 

44 

Trabzon 

39.43 

41.00 

30 

Table  2 

Recommended  average  day  and  declination  for  each  month  [31]. 

Months 

Day  of  the  year 

Date 

Declination  («5) 

January 

17 

17  Jan. 

-20.92 

February 

47 

16  Feb. 

-13.29 

March 

75 

16  Mar 

-2.42 

April 

105 

15  Apr. 

9.41 

May 

135 

15  May 

18.79 

June 

162 

11  June 

23.09 

July 

198 

17  July 

21.18 

August 

228 

1 6  Aug. 

13.45 

September 

258 

15  Sept. 

2.22 

October 

288 

15  Oct. 

-9.60 

November 

318 

14  Nov. 

-18.91 

December 

344 

10  Dec 

-23.05 
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Table  3 

Monthly  average  daily  global  (measured),  diffuse  (calculated)  and  extraterrestrial  solar  radiation  on  a  horizontal  surface,  clearness  index,  optimum  tilt-angle  and  monthly 
average  daily  global  radiation  on  optimum  tilted  surface  in  eight  provinces  of  Turkey  (radiations  are  given  as  MJm_2-day). 


Location 


Adana 


Ankara 


Diyarbakir 


Erzurum 


Istanbul 


Months 

H 

Hi 

H0 

I<T 

S„p,  (  ) 

Ht 

Jan. 

7.78 

3.79 

17.12 

0.45 

58 

12.43 

Feb. 

10.82 

4.87 

22.23 

0.49 

50 

14.80 

Mar. 

14.58 

6.20 

28.65 

0.51 

36 

16.93 

Apr. 

17.74 

7.62 

35.12 

0.51 

19 

18.38 

May 

21.20 

8.34 

39.50 

0.54 

4 

21.24 

Jun. 

23.16 

8.45 

41.22 

0.56 

0 

23.16 

jul. 

22.73 

8.24 

40.28 

0.56 

0 

22.73 

Aug. 

20.51 

7.57 

36.74 

0.56 

13 

20.89 

Sep. 

17.37 

6.35 

30.94 

0.56 

30 

19.36 

Oct. 

13.09 

5.05 

24.07 

0.54 

47 

17.42 

Nov. 

9.01 

4.00 

18.30 

0.49 

57 

14.20 

Dec. 

6.69 

3.47 

15.70 

0.43 

60 

10.97 

Jan. 

6.16 

3.36 

15.32 

0.40 

59 

9.74 

Feb. 

9.28 

4.55 

20.57 

0.45 

51 

12.79 

Mar. 

13.37 

5.98 

27.35 

0.49 

38 

15.73 

Apr. 

16.61 

7.54 

34.39 

0.48 

20 

17.32 

May 

20.87 

8.35 

39.32 

0.53 

6 

20.95 

Jun. 

23.35 

8.44 

41.32 

0.57 

0 

23.35 

Jul. 

23.88 

7.87 

40.26 

0.59 

2 

23.88 

Aug. 

21.25 

7.17 

36.24 

0.59 

16 

21.83 

Sep. 

17.45 

5.92 

29.85 

0.58 

33 

19.96 

Oct. 

11.86 

4.80 

22.53 

0.53 

49 

16.16 

Nov. 

7.36 

3.66 

16.53 

0.45 

58 

11.60 

Dec. 

5.11 

2.98 

13.88 

0.37 

60 

8.20 

Jan. 

7.93 

3.64 

16.54 

0.48 

59 

13.08 

Feb. 

11.49 

4.62 

21.70 

0.53 

51 

16.31 

Mar. 

15.56 

5.87 

28.24 

0.55 

38 

18.45 

Apr. 

19.81 

7.10 

34.90 

0.57 

20 

20.71 

May 

24.14 

7.45 

39.45 

0.61 

5 

24.20 

Jun. 

27.85 

6.61 

41.26 

0.67 

0 

27.85 

Jul. 

26.94 

6.58 

40.29 

0.67 

0 

26.94 

Aug. 

24.39 

6.02 

36.59 

0.67 

15 

25.03 

Sep. 

20.71 

4.87 

30.60 

0.68 

33 

23.92 

Oct. 

14.75 

4.33 

23.58 

0.63 

50 

20.77 

Nov. 

9.78 

3.68 

17.73 

0.55 

59 

16.40 

Dec. 

6.83 

3.34 

15.12 

0.45 

61 

11.60 

Jan. 

8.86 

3.07 

15.33 

0.58 

64 

17.08 

Feb. 

12.58 

3.89 

20.58 

0.61 

55 

19.74 

Mar. 

15.97 

5.44 

27.36 

0.58 

40 

19.62 

Apr. 

16.95 

7.51 

34.39 

0.49 

20 

17.70 

May 

19.89 

8.52 

39.32 

0.51 

6 

19.96 

Jun. 

23.17 

8.49 

41.32 

0.56 

0 

23.17 

Jul. 

23.19 

8.10 

40.26 

0.58 

2 

23.20 

Aug. 

21.21 

7.18 

36.24 

0.59 

16 

21.79 

Sep. 

17.25 

5.99 

29.86 

0.58 

33 

19.69 

Oct. 

12.57 

4.65 

22.54 

0.56 

50 

17.49 

Nov. 

8.79 

3.51 

16.54 

0.53 

61 

15.16 

Dec. 

7.03 

3.01 

13.89 

0.51 

65 

13.44 

Jan. 

5.00 

3.08 

14.70 

0.34 

58 

7.61 

Feb. 

7.75 

4.36 

20.00 

0.39 

50 

10.38 

Mar. 

11.82 

5.95 

26.89 

0.44 

37 

13.80 

Apr. 

15.76 

7.53 

34.12 

0.46 

21 

16.47 

May 

20.44 

8.41 

39.24 

0.52 

8 

20.54 

Jun. 

22.92 

8.56 

41.34 

0.55 

0 

22.92 

Jul. 

22.69 

8.23 

40.23 

0.56 

3 

22.71 

Aug. 

19.15 

7.65 

36.05 

0.53 

16 

19.68 

Sep. 

15.79 

6.23 

29.46 

0.54 

34 

18.01 

Oct. 

10.17 

4.86 

21.99 

0.46 

48 

13.57 

Nov. 

6.16 

3.47 

15.93 

0.39 

57 

9.45 

Dec. 

4.20 

2.70 

13.26 

0.32 

59 

6.55 

Jan. 

7.58 

3.55 

16.13 

0.47 

60 

12.58 

Feb. 

10.42 

4.66 

21.32 

0.49 

51 

14.50 

Mar. 

14.78 

5.94 

27.94 

0.53 

38 

17.48 

Apr. 

18.27 

7.41 

34.73 

0.53 

20 

19.06 

May 

22.39 

8.02 

39.41 

0.57 

6 

22.45 

Jun. 

25.31 

7.78 

41.28 

0.61 

0 

25.31 

Jul. 

24.90 

7.51 

40.28 

0.62 

1 

24.90 

Aug. 

22.35 

6.87 

36.47 

0.61 

15 

22.91 

Sep. 

18.64 

5.70 

30.35 

0.61 

33 

21.28 

Oct. 

13.35 

4.68 

23.22 

0.57 

49 

18.41 

Nov. 

8.83 

3.74 

17.32 

0.51 

59 

14.45 

Izmir 
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Table  3  ( continued ) 


Location 

Months 

H 

Hd 

Ho 

Kt 

$opt  (  ) 

Ht 

Dec. 

6.20 

3.24 

14.69 

0.42 

61 

10.34 

Samsun 

Jan. 

5.73 

3.14 

14.34 

0.40 

60 

9.39 

Feb. 

8.28 

4.34 

19.66 

0.42 

51 

11.41 

Mar. 

11.38 

5.88 

26.62 

0.43 

37 

13.20 

Apr. 

15.33 

7.51 

33.95 

0.45 

21 

16.01 

May 

19.44 

8.54 

39.19 

0.50 

7 

19.53 

Jun. 

22.59 

8.64 

41.35 

0.55 

0 

22.59 

jul. 

23.00 

8.14 

40.21 

0.57 

3 

23.02 

Aug. 

19.96 

7.43 

35.93 

0.56 

17 

20.56 

Sep. 

15.14 

6.28 

29.23 

0.52 

33 

17.17 

Oct. 

10.00 

4.79 

21.68 

0.46 

48 

13.31 

Nov. 

6.87 

3.45 

15.57 

0.44 

59 

11.19 

Dec. 

4.90 

2.80 

12.90 

0.38 

62 

8.34 

Trabzon 

Jan. 

5.29 

3.10 

14.45 

0.37 

59 

8.27 

Feb. 

7.81 

4.32 

19.76 

0.39 

50 

10.49 

Mar. 

10.52 

5.84 

26.70 

0.39 

35 

11.99 

Apr. 

13.62 

7.46 

34.00 

0.40 

19 

14.11 

May 

16.57 

8.66 

39.21 

0.42 

7 

16.63 

Jun. 

18.40 

9.15 

41.34 

0.45 

0 

18.40 

jul. 

16.71 

8.87 

40.22 

0.42 

3 

16.72 

Aug. 

14.44 

7.89 

35.97 

0.40 

14 

14.69 

Sep. 

12.32 

6.47 

29.30 

0.42 

30 

13.51 

Oct. 

8.84 

4.78 

21.77 

0.41 

46 

11.26 

Nov. 

6.12 

3.42 

15.67 

0.39 

57 

9.35 

Dec. 

4.56 

2.75 

13.01 

0.35 

61 

7.41 

Sopt  =  34.784-1. 3621(<5)-0.0081((5)2  (11) 

Sopt  =  34.783-1 .43 1 7(<5)-0.0081((S)2  +  0.0002(<5)3  (12) 

The  values  of  the  optimum  tilt  angle  calculated  from  Eqs. 
(10)-(12)  are  given  in  Table  4.  In  addition,  it  is  noticeable  that  the 
optimum  tilt  angle  for  the  month  of  June  is  negative;  the  negative 
sign  determines  the  orientation  of  the  solar  collector,  which  means 
that  the  solar  collector  is  faced  towards  the  north.  A  positive  sign 
indicates  that  the  solar  collector  is  directed  towards  the  south  [27], 

Table  5  shows  the  statistical  test  results  of  the  obtained  correlation 
equations.  It  can  be  observed  that  the  agreement  is  good  among  the 
results  obtained  by  Eqs.  (10)— (12).  The  optimum  tilt  angles  to  be  used 
in  winter  (December,  January  and  February),  in  spring  (March,  April, 
and  May),  in  summer  (June,  July  and  August)  and  in  autumn 
(September,  October  and  November)  for  each  province  of  Turkey 
are  given  in  Table  6.  Also,  the  values  of  latitude  +  15  and  the  yearly 
averages  are  presented  in  Table  6.  These  will  be  the  optimum  fixed 
tilt  angles  throughout  the  year.  The  values  of  the  yearly  total  global 
radiation  (in  MJ/m2-year)  on  a  tilted  surfaces  in  the  eight  provinces  of 
Turkey  are  given  in  Table  7  when  the  tilt  angles  of  the  solar  collector 
are  adjusted  by  taking  into  account  optimum,  seasonal,  yearly 
average,  latitude,  lat.+  15  and  lat.  — 15. 


5.  Statistical  methods 

In  the  literature,  numerous  statistical  methods  have  been  used 
to  compare  the  models  of  solar  radiation  estimation.  In  the 
present  study,  the  model  equations  are  evaluated  by  using  the 
most  widely  used  statistical  indicators  which  are  Mean  Bias  Error 
(MBE),  Root  Mean  Square  Error  (RMSE),  t-statistic  (t-stat)  and 
correlation  coefficient  (r)  [7-9,47-55]. 

5.1.  Mean  bias  error  (MBE) 

The  MBE  can  be  calculated  using  the  following  equation  [47] 

mbe=sE(K-*0  (13) 

i  =  1 


where  x,  is  the  i-th  measured  value,  yt  the  i-th  calculated  value, 
respectively  and  n  is  the  total  number  of  observations. 

The  MBE  provides  information  on  the  long-term  performance 
of  an  equation.  A  positive  value  gives  the  average  amount  of 
overestimation  in  the  estimated  values  and  vice  versa.  The  smaller 
the  value  of  MBE,  the  greater  the  performance  of  the  equation. 


5.2.  Root  mean  square  error  (RMSE) 


The  RMSE  can  be  computed  using  the  following  equation  [47] 


RMSE  = 


1/2 


(14) 


The  RMSE  provides  information  on  the  short-term  perfor¬ 
mance  of  an  equation.  The  value  of  the  RMSE  is  always  positive, 
representing  zero  in  the  ideal  case. 


5.3.  t-statistic  method  (t-stat) 


The  t-statistic  is  defined  as  follows  [48] 

stat=(  (n-PMBE2  V/2 
'vRMSE2-MBE2y 


(15) 


The  smaller  the  value  of  t-stat,  the  greater  the  performance  of  the 
model.  In  order  to  determine  whether  a  model’s  estimates  are 
statistically  significant,  one  simply  has  to  determine  a  t-critic  value 
obtainable  from  standard  statistical  tables,  i.e.,  ta/2  at  the  a.  level  of 
significance  and  (n- 1)  degrees-of-freedom.  For  the  model’s  estimates 
to  be  judged  statistically  significant  at  the  1  -  a  confidence  level,  the 
calculated  t-stat  value  must  be  less  than  the  t-critic  value  [48]. 


5.4.  Correlation  coefficient  (r) 

The  r  value  can  be  used  to  determine  the  linear  relationship 
between  the  measured  and  estimated  values,  which  can  be  calculated 
from  the  following  equation  [54] 

r  m  [(St— Sr)/St]1/2 


(16) 
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where  St  is  the  standard  deviation  and  Sr  is  the  deviation  of  the 
calculated  value  from  the  measured  value.  St  and  Sr  are  defined  as 
follows: 

n 

St=^(xa-x,)2  (17) 

i  =  1 

Sr  =  (x>-y>)2  <18> 


a 


Tilt  angle  (°) 


Tilt  angle  (°) 


Tilt  angle  (°) 


Tilt  angle  (°) 


where  xa  is  the  average  of  the  measured  values  and  is  given  by 

<19> 

In  this  study,  the  three  models  are  compared  on  the  basis  of 
the  above  mentioned  statistical  error  tests  and,  the  accuracy  of 
the  estimated  data  for  the  three  models  is  determined  using  these 


b 


Tilt  angle  (°) 


Tilt  angle  (°) 


Tilt  angle  (°) 


Fig.  1.  Monthly  average  daily  solar  radiation  availability  of  tilted  surfaces,  (a)  January-June;  (b)  July-December. 
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a  b 


tests.  For  better  data  modeling,  these  statistics  should  be  closer  to 
zero,  but  r  should  approach  to  1  as  closely  as  possible. 

6.  Mathematical  optimization  formulation 

The  maximization  problem  of  the  average  daily  radiation  on  a 
tilted  surface  (Hr)  can  be  expressed  as  an  optimization  problem  as 

maxHj- =  RH 


subject  to 

0<s<90  (20) 

where  HT  is  called  the  objective  function  to  be  maximized.  R  is  the 
functions  that  depend  on  the  design  variable  s,  while  H  is  a 
constant  for  a  given  month.  Since  the  optimization  problem  given 
in  Eq.  (20)  has  only  one  design  variable;  i.e.,  s,  the  optimum 
results  can  be  obtained  by  utilizing  either  the  graphical 
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Fig.  2.  Monthly  average  daily  solar  radiation  for  determined  tilt  angles. 


optimization  technique  or  one  of  the  numerical  optimization 
techniques  for  unconstrained  optimization  problems.  The  graphi¬ 
cal  representation  of  the  objective  function  is  given  in  Fig.  5. 

The  optimum  value  for  s  can  be  visually  read  from  the  figure. 
More  precise  results  for  s  can  be  obtained  using  one  of  the 
numerical  optimization  methods.  In  this  regard,  a  MATLAB  [56] 
based  program  has  been  written,  in  which  the  MATLAB  optimiza¬ 
tion  toolbox  was  used  to  solve  the  optimization  problem  given  in 
Eq.  (20).  An  unconstrained  optimization  command  in  the 
MATLAB,  fminbnd,  was  run,  which  is  based  on  golden  section 
search  and  parabolic  interpolation  [57],  The  optimum  value  of  s 


for  January  in  the  province  of  Adana,  for  example,  was  obtained  as 
58.32  and  the  corresponding  value  of  HT  was  12.43.  Consequently, 
it  can  be  seen  that  the  values  obtained  from  the  MATLAB 
optimization  toolbox  are  identical  to  the  values  found  by  the 
methodology  used  in  the  study  (Table  8). 


7.  Conclusions 

The  optimum  tilt  angles  of  the  solar  panels  may  be  determined 
when  the  tilt  angle  is  varied  from  0°  to  90  in  steps  of  1  °.  The  monthly 
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so 


Fig.  3.  General  correlations  developed  for  optimum  tilt  angles  of  solar  collectors 
used  in  Turkey,  (a)  Linear,  (b)  Quadratic  and  (c)  Third  order  polynomial. 

and  seasonal  optimum  tilt  angles  for  the  eight  provinces  having 
principal  climatic  conditions  of  Turkey  were  obtained  and  the  general 
correlation  equations  giving  the  monthly  optimum  tilt  angles  for 
Turkey  were  developed. 

According  to  the  above  results,  the  following  conclusions  can 
be  drawn: 

1.  It  was  found  that  the  optimum  tilt  angle  changed  between  0° 
(June)  and  65°  (December)  throughout  the  year  in  Turkey.  The 
optimum  tilt  angle  decreased  to  a  minimum  of  0  in  June  and 
July  for  Adana  and  Diyarbakir,  and  only  in  June  in  other 
provinces.  The  monthly  average  daily  total  radiation  at  this 
angle  was  generally  at  a  maximum. 

2.  The  optimum  tilt  angles  increased  during  the  winter  months 
and  reached  a  maximum  in  December  for  all  provinces.  The 


Fig.  4.  Comparison  of  the  results  obtained  in  the  literature  for  the  province 
of  Izmir. 


Table  4 

Values  of  the  optimum  tilt  angle  (in  degrees)  calculated  from  correlations 
obtained  for  Turkey. 


Months 

Eq.  (10) 

Eq.  (11) 

Eq.  (12) 

Jan. 

61.1 

59.7 

59.4 

Feb. 

50.7 

51.5 

51.9 

Mar. 

35.9 

38.0 

38.2 

Apr. 

19.8 

21.2 

20.8 

May 

7.0 

6.3 

6.3 

Jun. 

1.1 

-1.0 

-0.1 

jul. 

3.8 

2.3 

2.7 

Aug. 

14.3 

15.0 

14.5 

Sep. 

29.6 

31.7 

31.6 

Oct. 

45.7 

47.1 

47.6 

Nov. 

58.4 

57.6 

57.6 

Dec. 

64.0 

61.9 

61.0 

Table  5 

Statistical  test  results  (t-critic=3.106  for  tx=0.01  and  fc=12)  according  to  the 
correlation  equations  obtained  in  general  of  Turkey. 


Location 

Eq.  no 

MBE 

RMSE 

t-statistic 

r 

Adana 

(10) 

1.458 

2.182 

2.978 

0.9952 

(11) 

1.456 

1.754 

4.937 

0.9969 

(12) 

1.461 

1.659 

6.154 

0.9972 

Ankara 

(10) 

-0.042 

2.165 

0.064 

0.9952 

(ID 

-0.044 

1.054 

0.139 

0.9989 

(12) 

-0.039 

0.893 

0.147 

0.9992 

Diyarbakir 

(10) 

0.041 

2.391 

0.057 

0.9945 

(ID 

0.039 

1.382 

0.094 

0.9982 

(12) 

0.044 

1.316 

0.111 

0.9983 

Erzurum 

(10) 

-1.709 

2.728 

2.665 

0.9934 

(11) 

-1.711 

2.404 

3.359 

0.9949 

(12) 

-1.706 

2.462 

3.187 

0.9946 

Istanbul 

(10) 

0.041 

2.248 

0.208 

0.9935 

(11) 

0.039 

1.477 

0.088 

0.9976 

(12) 

0.044 

1.369 

0.106 

0.9979 

Izmir 

(10) 

-0.125 

2.004 

0.208 

0.9961 

(11) 

-0.127 

1.017 

0.418 

0.9990 

(12) 

-0.123 

0.963 

0.427 

0.9991 

Samsun 

(10) 

-0.542 

1.698 

1.117 

0.9971 

(11) 

-0.544 

0.983 

2.203 

0.9990 

(12) 

-0.539 

1.094 

1.881 

0.9988 

Trabzon 

(10) 

0.875 

1.288 

3.067 

0.9983 

(ID 

0.873 

1.449 

2.501 

0.9978 

(12) 

0.877 

1.394 

2.686 

0.9980 
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Table  6 

Latitude  +15,  seasonal  and  yearly  average  tilt  angles  for  provinces  of  Turkey. 


Location 

Lat. -15 

Lat.  +  15 

Spring 

Summer 

Autumn 

Winter 

Yearly  average 

Adana 

22.01 

52.01 

19.67 

4.33 

44.67 

56.00 

31.2 

Ankara 

24.56 

54.56 

21.33 

6.00 

46.67 

56.67 

32.7 

Diyarbakir 

22.55 

52.55 

21.00 

5.00 

47.33 

57.00 

32.6 

Erzurum 

24.55 

54.55 

22.00 

6.00 

48.00 

61.33 

34.3 

istanbul 

26.01 

56.01 

22.00 

6.33 

46.33 

55.67 

32.6 

Izmir 

23.25 

53.25 

21.33 

5.33 

47.00 

57.33 

32.8 

Samsun 

26.17 

56.17 

21.67 

6.67 

46.67 

57.67 

33.2 

Trabzon 

26.00 

56.00 

20.33 

5.67 

44.33 

56.67 

31.8 

Average 

24.39 

54.39 

21.17 

5.67 

46.38 

57.29 

32.56 

Table  7 

Yearly  total  global  radiation  (as  MJm~2-year)  on  a  tilted  surface  in  eight  provinces 
of  Turkey  according  to  the  adjusted  as  optimum,  seasonal,  yearly,  latitude,  lat.+15 
and  lat.  — 15  tilt  angles  of  the  solar  collector. 


Location 

Optimum 

Seasonal 

Yearly 

Latitude 

Lat.  +  15 

Lat. -15 

Adana 

6483.11 

6418.55 

6140.32 

6097.49 

5788.87 

6118.71 

Ankara 

6149.00 

6091.89 

5847.49 

5780.00 

5437.35 

5852.79 

Diyarbakir 

7481.56 

7401.19 

7052.27 

7006.02 

6635.73 

7027.77 

Erzurum 

6954.04 

6876.39 

6535.64 

6516.38 

6246.40 

6466.71 

istanbul 

5544.78 

5496.67 

5293.50 

5201.99 

4857.26 

5312.17 

Izmir 

6824.67 

6753.94 

6453.27 

6404.08 

6059.82 

6439.60 

Samsun 

5667.54 

5616.77 

5386.48 

5310.44 

4983.71 

5395.53 

Trabzon 

4661.38 

4623.18 

4450.55 

4386.25 

4131.21 

4453.94 

Table  8 

Values  of  Sopr  and  HT  of  each  province  for  January. 


Location 

Mathematical  optimization 

Method  used  in  the  study 

$opt  (  ) 

Ht 

$opt  (  ) 

Ht 

Adana 

58.32 

12.43 

58 

12.43 

Ankara 

58.76 

9.74 

59 

9.74 

Diyarbakir 

59.45 

13.08 

59 

13.08 

Erzurum 

63.94 

17.08 

64 

17.08 

istanbul 

57.82 

7.61 

58 

7.61 

izmir 

59.86 

12.58 

60 

12.58 

Samsun 

60.36 

9.39 

60 

9.39 

Trabzon 

58.76 

8.27 

59 

8.27 

monthly  solar  energy  collected  in  December  was  10.97,  8.20, 
11.60, 13.44,  6.55, 10.34,  8.34  and  7.41  MJm~2-day  for  Adana, 
Ankara,  Diyarbakir,  Erzurum,  Istanbul,  Izmir,  Samsun  and 


Trabzon,  respectively.  However,  these  values  of  the  monthly 
solar  energy  were  the  smallest  of  those  obtained  during 
the  year. 

3.  When  the  monthly  optimum  tilt  angle  was  used,  (i)  the  largest 
amount  of  monthly  solar  energy  collected  during  the  year 
were  23.16  (June),  23.88  (July),  27.85  (June),  23.20  (July), 
22.92  (June),  25.31  (June),  23.02  (July)  and  18.40  (June)  in 
MJm~2-day,  (ii)  the  yearly  collected  solar  energy  were 
6483.11,  6149.00,  7481.56,  6954.04,  5554.78,  6824.67, 
5667.54  and  4661.38  in  MJ/m2-year  for  Adana,  Ankara,  Diyar¬ 
bakir,  Erzurum,  Istanbul,  Izmir,  Samsun  and  Trabzon,  respec¬ 
tively.  This  value  (6824.67  MJm~2-year)  was  calculated  for 
Izmir  above  and  found  to  be  6820.36  MJm~2-year  by  Ref.  [18]. 

4.  When  the  monthly  optimum  tilt  angles  were  used,  the  yearly 
amount  of  collected  solar  energy  ranged  from  5635.86  to 
6483.11  in  MJ/m2-year  for  Adana.  Therefore,  the  enhance¬ 
ment  of  the  yearly  collected  solar  energy  was  1 5.03%.  These 
percentage  rates  were  14.07,  16.63,  21.58,  12.22, 15.83,  14.16 
and  13.00  for  Ankara,  Diyarbakir,  Erzurum,  Istanbul,  Izmir, 
Samsun  and  Trabzon,  respectively. 

5.  The  yearly  average  optimum  tilt  angles  for  a  south-facing 
solar  collector  were  31.2°,  32.7°,  32.6°,  34.3°,  32.6°,  32.8  and 
33.2°  in  Adana,  Ankara,  Diyarbakir,  Erzurum,  Istanbul,  Izmir, 
Samsun  and  Trabzon,  respectively.  When  the  yearly  average 
tilt  angles  were  used,  the  yearly  collected  solar  energy  were 
6140.32,  5847.49,  7052.27,  6535.64,  5293.50,  6453.27, 
5386.48  and  4450.55  in  MJm~2-year  for  Adana,  Ankara, 
Diyarbakir,  Erzurum,  Istanbul,  Izmir,  Samsun  and  Trabzon, 
respectively. 

6.  When  the  seasonal  tilt  angles  were  used,  the  yearly  collected 
solar  energy  was  6418.55,  6091.89,  7401.19,  6876.39, 
5496.67,  6753.94,  5616.77  and  4623.18  in  MJiTr2-year  for 
Adana,  Ankara,  Diyarbakir,  Erzurum,  Istanbul,  Izmir,  Samsun 
and  Trabzon,  respectively.  The  optimum  tilt  angles  for  spring, 
summer,  autumn  and  winter  on  average  for  Turkey  were 
21.17°,  5.67°,  46.48°  and  57.29°,  respectively. 

7.  The  seasonal  average  tilt  was  calculated  by  finding  the 
average  value  of  the  tilt  angle  for  each  season.  The  imple¬ 
mentation  of  this  requires  the  collector  tilt  to  be  changed  four 
times  a  year.  To  increase  the  efficiency  of  solar  collectors,  the 
solar  collector  should  be  mounted  at  the  monthly  average 
tilt  angle. 

8.  It  is  noted  that  the  amount  of  collected  solar  energy  for  all 
provinces  was  generally  in  order  of  optimum,  seasonally  and 
yearly  average  (in  Table  7). 

9.  The  model  Eqs.  (10)-(12)  obtained  in  the  study  can  be  used  to 
estimate  the  optimum  tilt  angle  of  a  solar  collector  in  Turkey. 
According  to  the  results  of  the  statistical  tests  of  the  correla¬ 
tion  coefficient  (r),  Eq.  (12)  generally  gave  the  minimum 
errors  when  compared  with  Eqs.  (10)  and  (11).  Also,  the 
equation  having  the  best  determination  coefficients 
(R2  =  0.9958)  among  the  three  models  is  Eq.  (12). 
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10.  The  model  equations  giving  the  greatest  results  for  each 
province  differ  according  to  the  statistical  methods  (Table  5). 
Considering  the  MBE;  the  model  equations  giving  the  best 
results  are  Eqs.  (11)  and  (12)  for  the  provinces  of  Adana, 
Ankara,  Diyarbakir,  Erzurum,  Istanbul,  Izmir,  Samsun  and 
Trabzon,  respectively.  Considering  the  RMSE;  the  model  equa¬ 
tions  giving  the  best  results  are  the  Eqs.  (10)-(12)  for  the 
provinces  of  Adana,  Ankara,  Diyarbakir,  Erzurum,  Istanbul, 
Izmir,  Samsun  and  Trabzon,  respectively.  Considering  the 
t-statistic;  all  of  the  model  equations  are  statistically  signifi¬ 
cant  excluding  Eqs.  (11)  and  (12)  for  both  Adana  and  Erzurum. 
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